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We have derived new boundary conditions on wave function at tlie normal metal / superconductor 
(NS) interface beyond effective mass approximation. These conditions are based on tight-binding 
approach and enable one to formulate quantitative model for tunneling spectroscopy of supercon- 
ductors with complex non-parabolic energy spectra. The model is applied to superconductors with 
unconventional pairing and with multiband electronic structure. In the case of single band un- 
conventional superconductors this model provides known conductance formula (Phys. Rev. Lett. 
74 3451 1995), but with generalized definition of the normal-state conductance. Based on new 
boundary conditions, we have calculated conductance in normal metal / superconducting pnictide 
junctions for different orientations of the NS interface with respect to the crystallographic axes of 
the pnictides. The present approach provides the basis for quantitative tunneling spectroscopy in 
multi-orbital superconductors. 

PACS numbers: 74.20.Rp,74.70.Xa,74.45.+c,74.50.+r,74.55.-|-v 



I. INTRODUCTION 



Tunneling spectroscopy was extensively applied up to 
now to reveal important features of electronic properties 
of superconductors ji It was predicted long time ago^*^ 
that in a tunnel junction between a normal metal and a 
metallic BCS superconductor (NS) the differential con- 
ductance dl/dV is expressed by the bulk density of 
states. Later, Blonder, Tinkham and Klapwijk (BTK) 
formulated the model for dl/dV in NS junctions with 
arbitrary barrier transparency by solving Bogoliubov- 
de Gennes (BdG) equation and explicit calculation of 
the Andreev and normal reflections coefficients It was 
shown that in the regime of low transparency of NS in- 
terface the resulting dl/dV corresponds to the energy 
spectrum of local density of states in a bulk supercon- 
ductor. On the other hand, dl/dV in high transparency 
limit is controlled by the Andreev reflection. The origi- 
nal formulation of the BTK theory has provided the ba- 
sis for tunneling spectroscopy and has been widely ap- 
plied for many conventional superconducting junctions 
where the symmetry of superconductor is conventional 
spin-singlet s-wave. The magnitude of energy gap for 
conventional s-wave superconductor has been determined 
with high accuracy. However, the effects of anisotropy 
of pair potential in d-wave pairing were not included in 
the original BTK formula. Therefore, after the discov- 
ery of high-Tc cuprates, the extension of BTK formula 
become really needed. Such an extension to unconven- 
tional superconductors has been formulated in Refs. |^ 
and 0. It was shown that tunneling conductance of un- 



conventional superconducting junctions is not always ex- 
pressed by the bulk density of states due to the presence 
of surface Andreev bound states (ABS)i^ii It has been 
revealed by a theory of tunneling spectroscopy of spin- 
singlet d-wave superconductor-^ that zero bias conduc- 
tance peak (ZBCP) stems from the surface zero energy 
ABS frequently observed in the experiments of high-Tc 
cuprates 4 This ABS has a fiat band dispersion and orig- 
inates from the sign change of the pair potential on the 
Fermi surface. Due to this flat dispersion, ZBCP ubiq- 
uitously emerges in actual experimentsi^'i'^ A number 
of anomalous quantum phenomena like non-monotonic 
temperature dependence of Josephson current in high Tc 
cuprate stem from this ABS>i2rJ^ Further, theory of tun- 
neling spectroscopy of normal metal /spin-triplet p-wave 
superconductor junctions has been formulatedi^^— stim- 
ulating by the discovery of Sr2Ru04ii^ It was shown that 
the line shape of a tunneling conductance in chiral p-wave 
superconducting junctions has a broad ZBCP due to the 
ABS with linear dispersion.—"— 

At present, to clarify the pairing mechanism of iron- 
based pnictides is one of hot topics in condensed matter 
physics. Just after the discovery of superconductivity 
in pnictides, s± symmetry of the pair potential (order 
parameter) has been proposed where pair potential 
changes sign between electron and hole Fermi pockets. 
The glue of this pairing is provided by spin- fluctuations, 
which typically appear in strongly correlated supercon- 
ductors. On the other hand, s++-wave pairing symme- 
try due to orbital fluctuations has been proposed as an- 
other candidate of superconducting pairingi^ The pos- 
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sibility of the inter-orbital pairing was also discussedi^ 
Since iron pnictides are multi-orbital systems with mul- 
tiple orbital systems, a theory of tunneling spectroscopy 
applicable to multi-orbital superconductors is strongly 
needed. Besides pnictides, there are many new uncon- 
ventional superconductors with multiple Fermi surfaces 
such as doped topological insulator C\ixBi2Se3^^— 

However, formulation of microscopic theory of tunnel- 
ing spectroscopy in multi-band superconductors is highly 
nontrivial task. The most crucial point is the bound- 
ary conditions on wave functions at the NS interface. 
Araujo and Sacramento have presented a new way to de- 
scribe boundary conditions between single-band normal 
metal and multi-band superconducting systems using 
phenomenological approach based on analogy between 
quantum waveguide theory and interband scattering;^ 
This idea was applied to actual pnictide junctionsr^i^i 
but the basis of this theory is not fully microscopic since 
interband and intervalley scattering effect is not fully 
taken into account. Other theories devoted to the study 
of coherent transport in junctions of iron pnictides are 
also phenomenologicali^"— The formulation of bound- 
ary conditions on the wave function in multi-band sys- 
tems has not become clear up to now. 

In this paper, based on equations of tight-binding 
model, we obtain the boundary conditions on the wave 
functions for the contact between a normal metal and 
a multiband superconductor. Up to now, tight-binding 
approach has been used for the study of ABS in various 
superconductors 42r— Here, we obtain boundary condi- 
tions beyond the effective mass approximation in order to 
take into account complex nonparabolic and anisotropic 
spectrum of energy band in the normal state and un- 
conventional pairing in multi-band superconducting sys- 
tems. The obtained boundary conditions provide an ex- 
tension of tight-binding approach by Zhu and Kroemer— 
to a superconducting case. The approach^ is physically 
transparent, since the only assumption is the prolonga- 
tion of solutions of tight-banding model on one addi- 
tional site to the left (right) sides of an interface. We 
apply the derived boundary conditions to the calculation 
of charge conductance between a normal metal and an 
iron pnictide superconductor for different misorientation 
angles between crystallographic axes of a pnictide super- 
conductor and the interface. Application of our theory 
to single band unconventional superconducting junctions 
allows one to reproduce preexisting formula of tunneling 
spectroscopy of unconventional superconductors, where 
the transparency of the junction in the normal state has 
a different momentum dependence. Brief account of some 
results of this paper is given in Ref. [10. 



II. ONE-DIMENSIONAL MODEL 

In order to understand the essence of new bound- 
ary conditions, we first consider onc-dimensional model 
of normal metal and spin-singlct s-wave superconductor 



junction. We use a model Hamiltonian H of the ID 
chain of atoms, whose electronic states are described in 
the tight-binding approximation, where Cooper pair is 
formed on the same site: 



H — Hn + Hs + Hi , 



(1) 



n<0,(T 



Hs = [t (aj^ „a<^,„+i -I- h.c.) - ^sa^a^,„] 



> 1 , cr 



H, 



7 {aifiO'<y-i + ^-c-) , 



(3) 



(4) 



with creation (annihilation) operator a\ „ {acr,n) of an 
electron with spin a on n site, and pair potential A. t' it) 
and /.i5 {ijln) are hopping and chemical potential in nor- 
mal metal (superconductor), respectively. Hjsi, Hs and 
Hi arc hamiltonian in the normal metal (N), in the su- 
perconductor (S) and at the interface, respectively. Eq. 
p]) is diagonalized by introducing the following canonical 
transformation: 



(5) 



which is a generalization of the Bogoliubov 
transformation^ to the case of a discrete lattice. 
In Eq. ([5]) al _^{a^^a) are operators of creation 
(annihilation) of quasiparticles satisfying Fermi anti- 
commutation relations, and are wave functions 
in BdG equation. The discrete version of BdG equations 



for wave functions u,. 



has the form: 



(6) 



where /i„ ~ iin (ms) foi' n < {n > 1), and i„ = t', 
7 and t for n < — 1, n = and n > 1, respectively. 
Basically, Eq. together with the corresponding self- 
consistent equations for the pair potential A„, provides 
the description for any spatially-inhomogeneous problem 
with an arbitrary set of hopping parameters between sites 
t„. However, this problem can be solved only numeri- 
cally. In order to formulate a relevant simplified model 
for an NS junction which allows analytical solution, we 
assume that there is no inhomogeneity of pair potential 
in a superconductor. That means, A„ = A for n > 1 and 
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n=2 n= 1 n=0 



n = l 



n = 2 



n = 3 




Using the boundary conditions Eqs. ([TO)) and (|lip and 
the wave functions Eq. ([7]), one can obtain 6, c given in 
Eq. 0: 



6 = 



(7i exp(iqZ) — exp(ifcZ) 
exp(ifc^) — (Ti exp(— iqZ) ' 



FIG. 1. Schematic illustration of one-dimensional model. Left 
region (red filled circles) corresponds to that of normal metal 
with hopping parameter t , right region (blue circles) corre- 
sponds to that of superconductor with hopping parameter t. 
The hopping parameter at NS boundary is given by 7. 



c = 7(l + 6)A, 



(12) 



with (Ti = tt' /j'^, and expression for the transparency aN 
at the interface: 



An = for n < 0. The structure under consideration is 
depicted in Fig. 1. 

First, we consider the conductance in a normal metal 
/ normal metal junction by setting A = in a super- 
conducting region. The electron with energy E{= e^) 
is injected from the left side and it is scattered at the 
interface. Then the wave functions u^^n for the left side 
^7i{= u^n) with n < and the right side 4'„(= u^^) 
with n > 1 are given by 



$„ = cyi'p{iqnl) + b cxp(— ign/), 
= c exp(ifcriZ), 



(7) 



where I is the lattice constant in a normal metal and 
a superconductor (for clarity, we consider them to be 
equal, but this restriction is not required4^) The first 
and the second term of represent an incident and a 
normal reflected waves, respectively. \l/„ corresponds to 
a transmitted wave. Here, q and k are determined by the 
equation 2t' cos{ql) = fj-jy + E and 2tcos{kl) = ^is + E 
with q,k > 0, respectively. The coefficients b and c are 
determined by the boundary conditions. These boundary 
conditions were proposed by Zhu and Kroemer4^ Their 
method is not limited by the assumption of parabolic 
single-particle excitation spectrum based on the elec- 
tronic effective-mass concept. In their idea, one can ob- 
tain the boundary conditions by the shift of the location 
of the boundary. If we shift the boundary to the right, 
we obtain the Shrodinger equation 



E^n 



-^Ar$o 



(8) 



where $1 is obtained by the natural extension of $„ (n < 
0) given in Eq. ([7]). The Shrodinger equation without the 
shift of the boundary is given by 



E^o = -^Ar$o + + 7*1- 



(9) 



By subtracting Eq. (jS]) from Eq. (|9]), wc obtain the 
boundary condition 



(10) 



Similarly, if we shift the boundary to the left, we obtain 
the boundary condition 



aN{k,q) = 1- I H = 



2 2(71 [cos[(g — k)l] — cos[{q + k)l]] 



1 + crf - 2(71 cos[{q + k)l] 



(13) 

The boundary conditions Eqs. (|10p and (|TT1) provide 
the conservation of probability flow J across the interface: 



2t' 2t 
Jn<-1 = — Im($;+i$„) = J„>i = _Im(1';+i*„). 

(14) 

After introducing finite difference derivatives in the fol- 
lowing form: ip[ = (*i - *o)/^'02 = ('^'i - *o)/^ the 
boundary conditions Eqs. PU)) and PT|) lead to usual 
boundary conditionSf^ obtained early in the continuum 
hmit. It is necessary to note that these boundary con- 
ditions, written in the form of finite differences, coincide 
with the mostly used Harrison's boundary conditions^ 
only for the case tri = 1. This feature of discrete bound- 
ary conditions Eqs. PU)) and PT|) were also mentioned 
in Ref. gl. 

By extending the method^ to the case of supercon- 
ducting junctions, one can obtain from Eq. (jG)) the 
following four boundary conditions at the interface be- 
tween normal metal and spin-singlet s-wave supercon- 
ductor junctions (see Fig. 1): 





= 7*1 




= 7*1 


7*0 


= i*o, 


7*0 





(15) 



where 5'„(<E>n) and *«(*&«) are the wave functions u^^n 
and Viy,n for electron and hole in S (N), respectively. 
These wave functions are given by 



' <f>n = exp(ign^) -f- b exp(— ign/), 
<&„ = a exp{iqnl), 

= c Mexp(ifcn/) + d v exp{—iknl) , 
^„ = c V exp{iknl) -\- d u exp{-~iknl) . 



(16) 



7$o = i*o- 



(11) 



The wave functions of a normal metal and a supercon- 
ductor contain four unknowns a, b, c, d describing the 
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Andreev and normal reflected waves in a normal metal, 
and two transmitted waves in the superconductor (c and 
d), where c (d) corresponds to transmission process by 
electron-like (hole-like) quasiparticles. These four un- 
knowns (a, b, c, d) are uniquely defined by four boundary 

conditions Eq. (jisp . In Eq. ([T6| . q, q {k, k) are wave vec- 
tors in normal metal (superconductor), corresponding to 
the energy E. 

Although q and q are real numbers, k and k become 
complex when | E |< |A| is satisfied. One can show that 
obtained wave functions provide the conservation of the 
probability flow by postulating boundary conditions Eq. 

. The expression for the probability flow on a discrete 
lattice (Fig. 1) follows from the BdG equations on the 
sites of the crystal lattice Eq. ^ : 



(17) 



It is necessary to note that the condition for the con- 
servation of probability flow at the interface between N 
and S, having the form of discrete sums (differences) in 
the crystal lattice (Eqs. (|14p . ((T7)) ). can be written in a 
quadratic form in terms of the probability amplitudes to 
be in states with wave vectors q,—q,q,k,—k multiplied 
on the group velocities in these states: 



dSn. 
dp 



dp 

+\d\' 



dss , 



dp 



den : 
dp 

Ip — — ■ 



dp 



\p=k 



(18) 



Eq. (lisp is similar to the corresponding expression in 
Ref. 11^ 

The above consideration of the tight-binding approx- 
imation of one-dimensional model of the NS junctions 
corresponds to equilibrium situation with zero voltage at 
the boundary V = 0. However, it can be generalized 
to the case of a finite voltage y 7^ on the microcon- 
striction of atomic sizes with a characteristic size much 
smaller than the elastic lei and inelastic Im characteristic 
mean free paths. In such pure microconstriction electron 
transport can be considered as a transport on the inde- 
pendent transverse modes. The current flowing through 
one mode is determined by the difference between the 
incoming f^{E) and outgoing f^{E) flows of electrons 
in the normal metali^ 



i{v) = m j{r{E)-r{E)}dE, 



(19) 



where f^{E) = /o(-E' — eV), fo{E) - the equilibrium 
Fermi distribution, rji = e/{'Kh), and 



riE) = AiE){l - n-E)) + B{E)r{E) 

+ {CiE)+D{E))foiE). (20) 



In Eq. (Uni) A{E), B{E), C{E) and D{E) are probabil- 
ities of the Andreev reflection, normal reflection, trans- 
mission as a electron-like quasiparticle and as a hole-like 
quasiparticle, respectively. The probabilities A, B, C, D 
in Eq. pO|) are calculated from the boundary conditions 
Eq. ([T5|) and the expressions for the probability flow Eqs. 
([itl) . P7|) . In calculating the probabilities A,B,C,D the 
incoming quasiparticle states must be normalized so that 
the probability flow in these states, described by the Eqs. 
(El) I ([T7)) . is equal to unity. This normalization provides 
a thermodynamic equilibrium in the absence of voltage 
V = on the NS junction. 

For the majority of superconductors, the magnitudes 
of A and E are much smaller than those of t and t' and 
the following conditions 



|AA|«l,|A/i'|«l 



(21) 



are satisfled. This is the so-called quasiclassical approxi- 
mation. Then, the relations q ~ g ~ go and fc ~ fc ~ fcg 
are satisfied, where fco and qo are momenta at the Fermi 
surface satisfying 2t'cos((7oO ~ f^N and 2tcos(fco0 — Ms- 
The resulting amplitudes a and b arc given by 



2gir(cos[(go - ko)l] - cos[(go + ko)l]) 
A 

A 



(22) 



with 

r = A/{E + y/E^ - A2), cxp{iqol) = 5, exp(ifcoO = S 
and 

A = -(1 - aiS~S){l - aii)[l - (1 - aN{ko,qo)T% 
00 



where aN{ko,qo) is defined by Eq. (|13p 

Within this approximation, we can reproduce the BTK 
result^ 



HV) = m I UoiE - eV) - fo{E)}aiE)dE, (23) 
where 



aiE) = 1- 



aN[l + TN I F p +{aN-l) I r 
I 1 - (1 - ajv)F2 |2 



(24) 



This is the well-known formula^^ with extended definition 
of the transparency crjv (see Eq. at the N/S inter- 

face. 

In this section, we have studied one-dimensional model 
with spin-singlet s-wave superconductor as the simplest 
case. However, the case of the contact between a nor- 
mal metal and a superconductor with anisotropic sign- 
changing pair potential on the Fermi surface is of most 
interest. The next two sections will be devoted to the 
consideration of this situation. 
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III. TWO-DIMENSIONAL MODEL FOR THE 

CONTACT OF A NORMAL METAL AND 
D-WAVE SINGLE BAND SUPERCONDUCTOR 

In this section, we extend our approach to unconven- 
tional superconductor. We show one typical example of 
two-dimensional lattice model of unconventional super- 
conductor. BdG equations on sites of the lattice of the 
d-wave superconductor in the x — y plane have the fol- 
lowing form: 



n-\-l,7n 

+ hi^runi+l + *n.m-l) - MS*n,rn 
+ Ao(4'ri+l,m + ^n-l.m ^ ^'n.m+l — ^n.m-l) 



(25) 



n,7n 

— Aoi'^n+l.m + ^n-l,m — ^'ri,m-|-l — ^'n,m-l) 



where ii,t2 arc hopping amplitudes between orbitals 
on sites, n, m are the numbers of site in x- and y- 
direction, respectively. The value of Aq is the ampli- 
tude of the anisotropic pair potential corresponding to 
the considered d-wave superconducting pairing: A{k) ~ 
2Ao(cos kx cos /cy ) , where kx and ky is quasimomcntum 
perpendicular and parallel to the interface, respectively. 

The boundary conditions for the contact of a normal 
metal and d-wave superconductor, described by the Eq. 
(PS)) . in the quasiclassical limit Eq. (PT|) are the same 
as boundary conditions Eq. ([T5|) . For the case under 
consideration, in boundary conditions Eq. 

means the wave function of layer n of atoms of d-wavc 
superconductor in the x ~ y plane. Due to the transla- 
tional symmetry in the y-direction in the electron (hole) 
wave functions ^'„,m (^n,m) we can omit second sub- 
script (m) corresponding to the coordinate of an atom 
in a direction parallel to the boundary. It should be re- 
marked that within quasiclassical approximation, these 
boundary conditions are satisfied for any type of uncon- 
ventional superconductors. 

Let's consider the situation, when the misorientation 
angle between interface and crystallographic axes of su- 
perconductor is equal to 7r/4. In this case the current 
through the two-dimensional pure microconstriction be- 
tween a normal metal and d-wave superconductor is de- 
termined by the integration over the transverse quasimo- 
mcntum ky of Eq. (|23p : Ip{V) = ^72/ dkyI{V,ky), where 
772 = S/27r, S - characteristic size of microconstriction, 
with the following definition of cr{E): 



E-2 - Al) with A± ^ A{±k^,ky). 

The present result is nothing but the formula of un- 
conventional tunneling conductance by Tanaka and 
Kashiway a^'^'^^ with generalized definition of ctat, pre- 
sented by Eq. ([T^ . The surface ABS is generated when 
the denominator of (t{E) becomes zero for = 0. This 
condition is given by 



1 = rr. 

Here, we consider dxy-'wave pairing with A^ 



A_ 



1 - (1 - aAr)rr |2 



4 Ao sin fca; sin /cj, . After simple manipulation, the energy 
level of surface ABS becomes E = for any ky on the 
Fermi surface. The dispersionless surface ABS is gener- 
ated in this case. The topological origin of this fiat band 
ABS has been clarified recently! ^^1^^ 



IV. TWO-DIMENSIONAL MODEL FOR THE 
CONTACT OF A NORMAL METAL AND A TWO 
BAND SUPERCONDUCTING PNICTIDE 

A. Zero misorientation angle 

Consider the application of this method for the case of 
two-dimensional electron transport through the bound- 
ary of a normal metal and a superconducting pnictide. In 
superconducting pnictidcs, there are two kinds of Fermi 
surfaces. One is the hole- like Fermi surface around F- 
point, and the other is the electron-like Fermi surface 
around the zone boundary. The minimum model to re- 
produce these Fermi surfaces is two-band model consid- 
ering the dzx and dyz orbitals in iron Sd-orbitalsi^ In this 
model, there are four kinds of hopping parameters ti, t2, 
ts and ti. As shown in Fig. [21 ti (i2) is the intra-orbital 
hopping between dzx (dyz )-orbitals in nearest neighbor 
site, and i4 are intra-orbital and inter-orbital hopping 
between next nearest neighbor sites, respectively. For 
the pair potential, s± model with A = 4Ao cos kx cos ky 
and S++ model with A = 2Ao(cosA:a; -I- cosky) + Ai are 
proposedi^i^ These pair potentials correspond to intra- 
orbital pairing and do not depend on the type of orbital. 

Let us first study the case of zero misorientation angle 
of the crystallographic axes of a pnictide with respect to 
the interface as shown in Fig. 2. The hopping perpen- 
dicular and oblique hopping between the sites on the left 
side and dzx (dyz)-orbitals on the right side are described 
by 71 (72) and 7^,7", (72; 72); respectively. For simpfic- 
ity, we assume that the periods of the crystal lattices in a 
normal metal and a pnictide arc the same. In the follow- 
ing calculations, wc drop oblique hopping for simplicity. 
The mathematical formulation of the problem and solu- 
tions of BdG equations in the two-dimensional case are 
given in Appendix. 

Proceeding as well as in the derivation of boundary 
conditions in the ID-model Eq. (|15p . but taking into 



with r = A+/{E + JE^ - Al) and F = A-/{E 



account independent hopping on dxz and dy^-orbitals of 
a pnictide, we obtain the following boundary conditions 
for NS junction with zero misorientation angle: 
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m = l 



»i = 




m = -l 



n= I n=0 n=l n=2 



n = 3 



FIG. 2. Boundary of normal metal / superconducting pnic- 
tide junction without misorientation. Left region (orange cir- 
cles) corresponds to the region of normal metal with hopping 
parameter t'i,t'2, right region (sites with two d-orbitals) cor- 
responds to the region of a superconducting pnictide with 
hopping parameters t\,t2,t3,tA- 71 (72) and 7( and 7I' (72 
and 72 ) are perpendicular and oblique hopping across the 
boundary to d^^^^dy z)-orhit&\s of a pnictide, respectively. 



(27) 



i'l^i =7i*?+72«'f, 
i;l>i =7i^?+72,*?, 
7i$o (ti + 2t3 cos fcy)^'o + 2it4 sin fcy*(^ 
+ 2AoC(fc,)^?, 

7i$o = {h + 2^3 cos fcy)§o + 2*^4 sin ky^'^ 

72*0 = {t2 + 2tsCos ky)^^ + 2it4suiky^'^ 
+ 2AoC(fc,)*(^, 

72$o = {t2 + 2t3 COS ky)^^ + 2iti sin ky^'^ 
- 2AoC(fc,)vI'(^, 



where C,{ky) ~ cosky and 1/2 for s± and 5++ models, 
respectively. Due to the translational symmetry of this 
structure in a direction parallel to the boundary, ky com- 
ponent of the quasimomentum is conserved. Also because 
of this translational symmetry in the electron (hole) wave 

functions ^'n.m'' (^n.m''), thc sccoud subscript (m) corre- 
sponding to the coordinate of an atom in a direction par- 
allel to the boundary is omitted. The wave functions 
of thc NS contact arc defined by 6 plane waves with 
amplitudes a,b,ci,C2,di,d2- a, b describe the Andreev 
and normal reflected waves in normal metal, ci (02) and 
di (0^2) describe electron-like and hole-like transmitted 
waves corresponding to inner (outer) Fermi surface in a 
superconducting pnictide, respectively: 



= exp(iqinl) + b exp(— igin/), 

= a exp(ig2"-07 
4"" = ciMi(fci) exp(ifcin/) -t- 021*1(^2) exp{ik2nl) 
+ diuilks) exp(ifc3n/) -f- ^2^1(^4) exp(ifc4n/), 

= ciM2(fci) exp(ifcin/) + 02^2(^2) exp(ifc2"-0 
+ diU2{k3) cxjplik^nl) + d2U2{k4) exp(ifc4n/), 

= ciwi(fci) exp{ikinl) + C2Wi(fc2) exp{ik2nl) 
+ divi{k3) exp{ik3nl) + ^2^1(^4) exp(ifc4n/), 

^'fj = ciV2{ki) exp{ikinl) -I- 02^2(^2) exp(ifc2'^0 
+ diV2{k3) exp{ik3nl) + d2V2{k4) exp{ik4nl). 



(28) 





0.06 
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" 0.02 
0.01 
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FIG. 3. The excitation spectrum of a superconducting pnic- 
tide for a fixed value of ky. (A) - misorientation angle is 
equal to 0, ky = 0, (B) - misorientation angle is equal to 
7r/4, ky = 0, (C) corresponds to the (A), depicted on a larger 
scale, (D) corresponds to (B), depicted on a larger scale. Red 
dashed line and black solid line correspond to the two different 
bands. 
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FIG. 4. The Fermi surface of a pnictide in the unfolded Bril- 
louin zone. (A) - misorientation angle is equal to 0, (B) - 
misorientation angle is equal to 7r/4 

where qi, q2 are the wave vectors for electron and 
hole in normal metal with excitation energy E, respec- 
tively, fci (^2) and ks (^4) are the wave vectors for 
electron-like and hole-like quasiparticle corresponding to 



7 



inner (outer) Fermi surface in pnictides. Six coefficients 
a,b,Ci,C2,di,d2 in Eq. (pS)) are uniquely determined 
from six boundary conditions in Eq. (f27|) . Electron and 
hole coefficients Ui{kj) and Vi{k.j) in wave functions Eq. 
p8)) are also found from Eq. (jAip . The excitation spec- 
trum of a superconducting pnictide, corresponding to a 
fixed value of fcj, = at zero angle of misorientation is 
shown in Fig. 3 (A),(C). The corresponding Fermi sur- 
face is shown in Fig. |4l There are four intersection points 
with ky = 0. At these points, E(k^) has minima shown 
in Figs. 3 (A) and (C). The existence of the four quasi- 
particle states in a pnictide with certain sign of the group 
velocity follows from Figs. 3 (A) and (C). 

The expression for the probability flow with fixed wave 
vector ky in the direction parallel to the x axis follows 
from the BdG equations on sites of the crystal lattice of 
a pnictide Eq. (|A1[) and has the following form: 

Jp = ^{{h + 2t,cosky)lm{{n+im - 

+ (t2 + 2<3Cosfc,)Im{(*f+i)**^-(^'^+i)**^} 
+AUsmkyRe{{n+ir'^i + 

+2AoCOsfc,Im{(v[/^+i)*^^ 

+ K+iVH + (*f+i)**fJ)- (29) 

One can show that the boundary conditions Eq. ((27)) 
provide the conservation of the probability flow J ~ Jp 
across the interface between a normal metal and a super- 
conducting pnictide for each value of ky. As well as in the 
case of previously considered ID model (Eq. ((T8)) ). the 
condition of flow conservation at the NS boundary, hav- 
ing the form of discrete sums (differences) on sites of the 
crystal lattice (Eqs. p^ . (|29l) ). in the case of zero misori- 
entation angle between crystallographic axes of a pnictide 
and the interface can be written in a quadratic form of 
the amplitudes of the probability to be in states with 
quasi- momenta qi,q2,ki, i = 1..4 multiplied by group 
velocities in these states. 

The current through the two-dimensional microcon- 
striction between a normal metal and a supercon- 
ducting pnictide is determined by integration of Eq. 
(Uni) over the transverse quasimomentum ky-. Ip{V) = 
'72 / dkyI{V, ky). In this case the probabilities C, D of the 
quasiparticle propagation into a superconductor, given 
by Eqs. p9)) and ([20]), are determined by the sum of 
the scattering probabilities into individual bands: C ~ 
Ci + C2,D = D1+D2. The coefficients A, B,Ci,C2, Di, 
and D2 in Eq. (PO)) are calculated from the boundary 
conditions Eq. (|27|) and the expressions of the probabil- 
ity flow in Eqs. (fH)) and (P^ . In the actual calculations, 
one should take into account that the original quasipar- 
ticle states should be normalized so that the probability 
flow in these states, described by the Eqs. ([T4l). (p9|) . is 
equal to unity. 

It is possible to demonstrate, that taking into account 
oblique hopping between the boundary 71,72, 7", 72 (see 



Fig. 2) allows to obtain in quasiclassical limit Araiijo 
and Sacramento boundary conditions^ only for the spe- 
cial case, when the following relations between hopping 
amplitudes are fulfilled simultaneously: 




FIG. 5. 2D NS boundary. Angle between crystallographic 
axes of a pnictide and a normal metal is equal to 7r/4. The 
lower left region (orange circles) corresponds to the region 
of normal metal with hopping parameters t'i,t'2, right region 
(sites with two d-orbitals) corresponds to the region of a su- 
perconducting pnictide with hopping parameters t\,t2,t^,t4,. 
Boundary is described by hopping parameters 71 , 72 , 7i , 72 ■ 



B. Non-zero misorientation angle 

The proposed method allows one to consider the co- 
herent electron transport in NS structures with non-zero 
misorientation angle as well. It is necessary to note 
that the microscopic calculation of the conductance for a 
non-zero misorientation angle of a pnictide crystal with 
respect to the boundary is presented here for the first 
time. Previous phenomenological approaches^^*^ don't 
allow one to carry out such calculations. In considering 
the electron transport across NS contact with a nonzero 
misorientation angle, it is necessary to take into account 
hopping at the two adjacent atomic layers of a pnictide 
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(FiglS]). BdG equations in the considered case corre- 
sponding to s± symmetry of the pair potential in pnic- 
tides are given in the Appendix, sec Eq. ()A2p . Hop- 
ping across the NS boundary for non-zero misorienta- 
tion angle is described by a larger number of parameters, 
rather than at zero misorientation angle between crystal- 
lographic axes of a pnictide and the interface (see FiglS]). 
In addition to hopping parameters 71 and 72, we should 
use additional parameters of hopping across the bound- 
ary j[ and 72- These parameters of hopping across the 
boundary take into account connection of orbitals from 
the last atom layer of a pnictide with the penultimate 
from the boundary atom layer of the normal metal. Tak- 
ing into account these processes is necessary due to the 
breaking at the boundary of the diagonal bonds in the 
crystal lattice of a pnictide for non-zero angle of misori- 
entation (see Figl5]). Also in the normal metal together 
with the nearest neighbor hopping t[ we need to consider 
the diagonal hopping t'2 in square lattice. 

The wave functions Eq. (|A4[) and the relation for prob- 
ability flow (|A5|) take into account not only the electron 
transport in two energy bands, but also in two valleys 
in these bands (see Fig. [3j3,D). It is known from the 
physics of semiconductors, that interference of states in 
the valleys is possiblei^ This interference leads to the 
fact that the condition of flow conservation at the bound- 
ary of the NS contact, having the form of discrete sums 
(differences) on sites of the crystal lattice, in the case 
of nonzero misorientation angle between crystallographic 
axes of a pnictide and the interface can not be written 
in a quadratic form of the amplitudes of the probability 
to be in states with quasimomentum qi,q2,ki, i = 1..8, 
multiplied by group velocities in these states. 



C. Numerical results 

Here we will show the results of numerical calcula- 
tions of angle- resolved conductance(dI/dV) as a func- 
tion of bias voltage V in normal metal / superconduct- 
ing pnictide junctions. We use the following values of 
hopping parameters and chemical potential in a pnictide: 



(B) 
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0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 
eV 



FIG. 6. Angle resolved conductance without misorientation 
for (A) s±-model and (B) S-|-+-model. Value of the quasi- 
momentum, parallel to the interface, ky — 0.01. The values 
of hopping parameters at the interface are chosen as 71 — 
0.1,72 = 0.14 (eV) (dashed line), and 71 = 0.009,72 = 0.005 
(eV) (solid lines). 
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FIG. 7. Same as Fig. [Slbut with ky = Sn/A. 
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FIG. 8. Same as Fig. [6] but with ky = 5n/ 



h = -0.1051, t2 = 0.1472, ts = -0.1909, U = -0.0874 
and fis = -0.081 (eV), according to Ref. [H. We 
assume the momentum dependence of the pair poten- 
tial in the s± model as A±{k) ~ 4Ao cos fc^ cos fcj, with 
Ao = 0.008 (eV), and s++ model in the form A = 
2Ao(cos k^ + cos ky) + Ai with Aq = 0.002, Ai 0.0042 
(eV). In Figs. [HltolHl the magnitudes of tunneling conduc- 
tance normalized by their values in the normal state are 
shown for s± and s-|_+ models with zero misorientation 
angle. The hopping parameters and chemical potential 
in normal metal are t[ = 0.3, ti^ = 0, hn = 0.2. For 
the hopping parameters at the interface, we choose two 
cases with 71 = 0.009,72 = 0.005 (eV) (low transmis- 
sivity) and 71 = 0.1,72 = 0.14 (eV) (high transmissiv- 
ity). Calculated charge conductance dl/dV for low and 
high transparent junctions with ky = 0.01 is shown in 
Fig. El The horizontal axis represents eV normalized by 
Amaa; ; whcrc /S.max IS tlic maximum of two gaps for fixed 
ky. One can clearly sec two gap features reflecting the 
presence of two kinds of Fermi surfaces (see Fig. [3]) in 
both s±- FiglHUA) and s+-|-- cases Fig[6UB). In the case of 
ky ^ 0, the interorbital hopping is absent. Therefore, 
the obtained conductance can be represented by a simple 
summation of the individual orbital's contributions. On 
the other hand, in the case of s±-wave with low trans- 
missive interface, a sharp subgap peak appears in the 
energy gap as shown in Figl7]^A) and FiglHl^A), respec- 
tively. These energy structures do not correspond to the 
density of states in the bulk. Since dl/dV corresponds to 
the energy spectrum of local density of states in the low 
transmissivity, we can conclude that these subgap struc- 
tures originate from the surface Andreev bound states at 
finite energies. The bound states disappear in the case of 
high transparency of the interface. On the other hand, as 
is seen from FigjT^B) and FiglSKB), these features are not 
present in the case of s++-wave when the signs of pair 
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FIG. 9. Angle resolved conductance with misorientation 
angle 7r/4 for (A) s±-model and (B) s++-model. Value of 
the quasimomentum, parallel to the interface, ky — 0. The 
values of hopping parameters at the interface are chosen as 
71 = 0.1,72 = 0.14, 7( = 0.2,72 = 0.06 (eV) (dashed line), 
and 71 = 0.009, 72 = 0.005, 7^ = 0.02, 72 = 0.01 (eV) (solid 
lines) . 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 0.0 0.2 04 0.6 0.8 1.0 1.2 1.4 

eV eV 



FIG. 10. Same as Fig. [9] but with ky = n/3. 



potentials in different bands are the same. From these 
results, we can conclude that the surface Andreev bound 
states are formed in the s± case due to the sign change of 
pair potential and the interorbital hopping t^. Note that 
a sharp subgap peak in the angle-resolved conductance 
discussed here should be broadened after summation over 
ky is made, as calculated by Onari et al within different 
mo del 4^ 

Next, we calculate the case with finite misorientation 
angle 7r/4. In this case, we shall introduce additional 
hopping coefficients in a normal metal (^2) and at the 
interface (7^, 73) corresponding to the direction perpen- 
dicular to the interface. We choose t'2 = 0.01 (eV), 
7^ = 0.02,72 = 0.01 (eV) (low transmissivity case) and 
7^ = 0.2,72 = 0.06 (cV) (high transmissivity case). 
Other hopping parameters arc the same as in the case 
with zero misorientation. In Figs. |9] and [101 we show 
the angle resolved conductances for fcy = and 7r/3, re- 
spectively. In the case of s± model, one can see the two 
gap features for both ky = and 7r/3 as shown in Figs, 
inij A) and [TUI A) . Subgap peaks arc absent even in low 
transparent junctions, in contrast to the case with zero 
misorientation angle. This is because no sign change at 
fixed ky values occurs in the case when misorientation 
angle equals 7r/4 (Fig. HJ^B)). For the same reason, gap 
structure without subgap peaks appears also in the case 
of S++ model. 

Let us summarize the results of the conductance of 
normal metal / superconducting pnictide junctions. In 
the case of s++ model, only the two-gap structure with- 



out subgap peaks appears for any misorientation angle 
and any value of ky. On the other hand, in the case of 
s± model with low transparent interface, subgap peaks 
appear for zero misorientation angle and finite ky . These 
subgap peaks originate from sign change of the pair po- 
tential at fixed ky values in the presence of the interor- 
bital hopping. 



V. CONCLUSION 

In this paper, we have presented consistent tight- 
binding model for the coherent electronic transport in the 
contact between a normal metal and a superconductor. 
Based on a tight-binding model beyond effective mass 
approximation, we have derived boundary conditions on 
a wave function at a contact between a normal metal and 
a superconductor with unconventional pairing symmetry. 
We have extended the previous tight-binding approach 
used in semiconducting hcterostructures^ to the case 
of superconducting junctions. The obtained boundary 
conditions contain real space information only without 
any momentum derivatives, and they have clear physi- 
cal meaning. These conditions provide current conser- 
vation and enable one to formulate consistent approach 
for tunneling spectroscopy of superconductors with com- 
plex nonparabolic energy spectrum, including multiband 
electronic structure and unconventional symmetry of su- 
perconducting pairing. We have shown that application 
of this theory to single-band superconductor junctions 
allows one to reproduce the preexisting conductance for- 
mula 0. Based on the derived boundary conditions, we 
have calculated conductance in normal metal / supercon- 
ducting pnictide junctions for different misorientation an- 
gles between the interface and the crystallographic axes 
of a pnictide. The present approach provides the basis 
for tunneling spectroscopy of multi-orbital superconduc- 
tors. Moreover, this approach is suitable for a consistent 
description of electronic transport in structures with sur- 
face states described by Majorana fermions in topologi- 
cal superconductors , ^^1^°" — which would be the subject 
of our future study. 
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Appendix A: Derivation of equations of zero misorientation angle of the crystallographic axes of 

two-dimensional model pnictides with respect to the interface have the following 

form: 

Bogoliubov-de Gennes equations on sites of the Fe crys- 
tal lattice in the x — y plane of a pnictide for the case of 



Vl/° 



^si^n+l.m+l + 771-1 + ^n+l.Tji-l 

*n-l,r7i+l) + *4( — ^(^+i.„i+i — *n-l, 777.-1 



'n+l, 777-1 ^ ^77-1,777+1,' l^ti '^n,m "UV ^71 + 1,777 + 1 

' n— l.rn — 1 ' n+l^rn — l ' n — l,m+l/ ^ Ti,rn ' 
i2(*f+i,„ + + tl(*l7„+l + Km-l) 

+ ^3(^^+1, ,71+1 + "^n-l.m-l + *n+l,77i-l 
+ *n-l,77i+l) + *4(— ^"+l,,ri+l — '^n-l,m-l 
+ *n+l, 777-1 + *«-l, 777+1) - fJ-S'^i.m + Ao(*(^+i_„+i 



71—1.777—1 ' 77 + 1,777 — 1 ' 77 — 1,777+1 / ^^71,777' 



+1,777 



+ 1,777 + 1 



n— 1.771/ 



^ 77+1,771—1 



*7i-l,77i+l) + ^n+l,,ri+l '^n-l.m-l 

^n+l,m-l + ^n-1,777+1) ~ Ms*", 777 ~ ^o(*"+l, 777+1 



' ^77 + 1,777 — 1 ' ^ 77 — 1, 777+1 7 



=,/3 







V/3 



^3(^^+1, 7,1 + 1 T" ^77-1,777-1 ^ ^77+1,777-1 

■ ^n-l,7n+l) + ^4( — ^„+l,7n+l ~ ^77-1,777-1 

■ *n + l, 777-1 + ^«-l, 777+1) - A*S*^^,777 - ^0(*,''i + l,,77+l 



P 

n+l,m — 1 



p 



1) 



(Al) 



where i = 1..4, are hopping amplitudes between or- 
bitals on sites in a pnictide in the two-orbital model.— 
The value of Aq is the amplitude of the anisotropic 
pair potential corresponding to the considered s± super- 
conducting pairing model: A+(fc) = 4 Aq cos fc^; cos fcy 
fcy, kx are parallel and perpendicular to the interface com- 
ponents of quasimomentum respectively. The wave func- 
tions of a superconducting pnictide have the upper or- 
bital index q;(/3): ^'^^^\ corresponding to dxzidyz) or- 
bital respectively. The subscripts n, m of the wave func- 
tion "^nlm of a pnictide describe the coordinates of sites 
of the crystal lattice (Fig. 2). As well as in the con- 
sidered above ID-model ^'"fm in Eq. (jAlj) describe the 
electron states, and ^Z^S - hole states. 

For misorientation angle 7r/4 between crystallographic 
axes of a pnictide and the interface (Fig 15]), Bogoliubov- 
de Gennes equations on sites of the Fe crystal lattice in 



the X — y plane of a pnictide differ from Eq. (jAip and 
have the following form: 
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+ ^n-l,m+l) 
hi'^n+l.m+l + *n-l,m-l) 

hiK+2,n + K-2.n + K.rn-2 + *",™+2) 

-2 + ^n,m+2) 



n— 2,?n 

/3 



71 — 2 , rn 



^ *L„-2 
n,m — 2 



n.rnf 

r,/3 



Ms*f,.,„+i2(* 



-ii(*fi+i,„i+i 

" i4(-*"+2,m " 

r,/3 



n+1 ,m — 1 
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^— l,m+l 



) 



n— 2.??i 
n—2.ni 
71 — 2, m ' 



71.771 



+ *n,m+2) 
*ri,m-2 + *n,m+2) 

:) 



71,777 — 2 



^ 71,771' 

+ ^2(^'"+l,„i+l + *"_i,,„_i) 



^n,m+2 



,m 



n — 2.771 



n. — 2 , 777 



71 — 2,771 

-t4(-*l^+2.,„ - 

^ 71,771.' 

-Ms'J'l™+t2(*f+l,™ 

+ i2(*^+i,™+i + *^_i 

77 — 2,777 



n,7n — 2 

77,771 — 2 



77,777 + 2 J 



,,777 + 2/ 



+ ^71,777 - 2 + ^77,771 + 2) 

-1 + 777+1) 



,777— 1 > 



77,777- 

n.m — 2 



^n,m+2) 



+ ^3(^^+2,™ 
+ i4( — *"+2,m ~ *n-2.m 
- MK+2,m + K-2,m + K,m-2 + K,ra+2) 



a \ 
n, 771+2/ 



(A2) 



+ *^(7ie-*^-«' + 72e*^'''') + 7i*2 + 72*2 > 

+ *^(7ie-*'=«' + 726-*'="') + 7^*^ + 7^*^, 
$0(716*''''' + 72e"'''»') + 7i$-i = ti^'fje*''''' 

+ ^2*06"'''''' + i3*-i - t4*^!i + Ao^-^i, 

$0(716*''''' + 726"*''''') + 7l'l-i = ti^'EJe*''''' 

+ ^2*0 6"'''''' + is^-.l - - ^O'^-I, 

$0(716^*''''' + 72e'''»') + 72$-i = ti*o6*''''' 

$0(716^*''''' + 726*''''') + 72$-l = ^1*06*''''' 
+ t2^'o6-'^'''' + t3*-l - ti^-l - Ao*^i, 
7i*0 = ^3*0 - + Ao*^, 

7i$o = ^3*0 - i4*(^ - Ao*?, 

7^$o = i3*(^ - ^4*0 + Ao*^, 
7^$o = i3*(? - ^4*0 - Ao*^, 

i2^i=7i*?+72*f, 
[t^$i=7i*?+72^?- 

(A3) 



As in the previously considered case of boundary con- 
ditions for zero misorientation angle Eq. (|27|) . due to 
translational symmetry in the direction parallel to the 
boundary, in electron (hole) wave functions ^'n.m'' i^n[rn) 
second subscript (m) corresponding to the coordinate of 
an atom in a direction parallel to the boundary is omit- 
ted. 



The boundary conditions for the contact between a 
normal metal and a pnictide, considered in the frame- 
work of the two-orbital model, for misorientation angle 
7r/4 between crystallographic axes of a pnictide and the 
interface have the following form: 



The wave functions in a normal metal / supercon- 
ducting pnictide contact in the case of misorientation 
angle 7r/4 between crystallographic axes of a pnictidea 
and the interface are defined by eight plane waves with 
amplitudes ai, 61, 02, 62, ci, C2, di, (i2, /i, /2, ffi, 92- Here 
the coefficients ai, 61, 02, 62 describe Andreev and normal 
reflected waves, while ci, C2, di, ^2, /i, /2, ffi, ff2 describe 
eight waves transmitted into a two-band superconduct- 
ing pnictide: 
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' <&„ = exp(iqinl) + bi exp{—iqinl) + &2 exp(— ^(72?^0^ 

$„ = ai cxp(i(73nZ) +02 cxp(iq4riZ), 

= ciUi(fci) exp(ifcin/) + C2Mi(fc2) exp(ifc2nZ) 
+ diUi(fc3) expiik^nl) + d2Ui{k4) exp{ik4nl) 
+ fiui{kr^) cxp{ik5nl) + f2Ui{ke) cxp{ikenl) 
+ giui(k-j) eyi\i(ik'jnl) + g2Ui{ks) e:iq>{iksnl) , 

= ciU2(fci) exp(ifcin/) + C2M2(fc2) exp{ik2nl) 
+ diU2{k3) exp{ik3nl) + ^2^2(^4) exp(ikinl) 
+ fiU2{k5) exp{ik5Til) + /2W2(fc6) exp{ikGnl) 
+ ^1^2(^7) exp{ik7nl) + g2U2{ks) exp{iks'nl), 

4'" = ciwi(fci) exp(ifcinZ) + C2Wi(fc2) exp{ik2nl) 
+ divi (fca) exp(ifc3nZ) + ^2^1 (^4) exp^ik^nl) 
+ fiviih) expiiksnl) + f2Vi{ke) exp{ikenl) 
+ givi^kj) exp{ikTnl) + g2Vi{ks) exp{iksnl) , 

4'^ = CiV2{ki) exp{ikinl) + 02^2(^2) exp{ik2nl) 
+ diV2{k3) explik^nl) + d2V2{k4) expiik/^nl) 
+ /iW2(fc5) cxp{ikznl) + f2V2{kQ) cxp{ikenl) 
+ giV2{k-j) cxp{ik-jnl) + 52^2(^8) cxp{iksnl). 

(A4) 

Four transmitted waves with amplitudes ci, C2, di, d2 cor- 
respond to the lower band, depicted by black solid line on 
Fig. [3l[B),(D). These four waves are propagating waves 
except the energy range lower than the superconducting 
gap Aq. Four plane waves with amplitudes /i,/2,3i,52 
correspond to the upper band, depicted by red dashed 
line on Fig. |3l[B),(D). These four waves arc evanescent 



waves on the scale of pair potential Aq. 

Expression for the probability flow in the case of mis- 
orientation angle between crystallographic axes of a pnic- 
tide and the interface equal to 7r/4 differs from the cor- 
responding relation for the case of zero misoricntation 
angle Eq. ((29)) and has the following form: 

+t2lm{iK+irKe-^''''} 
+tslm{{n+im-i + (*"+2)**"} 
+hlm{{^t+irHe-"'^'} + t2lm{(vl'f+,)*vl,^e"=«'} 

-M(im{(vi'^+i)*vi/^_ J + iMK^irn-i} 

+Im{(vI/^+2)**^} + lm{{^^_,2m}) 
-tilm{(*^+,)**^e»'^«'} - t2lm{(i'5J+i)*i'«e-''=«'} 

-tiIm{(§f,+i)*§^e-^'=«'} - t2lm{{^^^^,r¥y^y^ 

+im{(§«+2)**^} + im{{^i^2rn}) 
+Aoim{(vi/^+i)**"-i + in+iTn-i 

+iK+2rK + i<+2r^i + (^^+2)**f })■ (A5) 
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